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DYNAMICS  OF  OBSERVED  REALITY:  ABRIDGED  VERSION  OF  CLASSICAL 
AND  QUANTUM  MECHANICS 


Azizul  Haque  tnd  Thomas  F.  George 

Departments  of  Chemistry  and  Physics  &  Astronomy 

239  Fronczak  Hall 

State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260  USA 


INTRODUCTION 

The  present  paper  is  concerned  with  a  unified  description  of 
classical  and  quantum  physics.  A  different  way  of  thinking  about  the 
quantum  world  is  expounded,  which  already  has  led  to  significant  results 
in  statistical  mechanics  and  holds  future  promise  for  quantum  mechanics. 
An  explanation  is  provided  for  the  original  motivation  in  quantum  theory, 
the  formalisms  that  have  evolved  from  it  and  their  differences  with  the 
classical  theories.  The  conceptual  difficulty  that  permeates  our  view  of 
the  microscopic  world  is  spelled  out  and  a  remedy  thereof  suggested. 


BACKGROUND 

The  laws  of  nature  for  the  microscopic  world  are  very  much  different 
from  those  of  the  macroscopic  world,  though  the  building  blocks  of  any 
macroscopic  system  are  the  microscopic  particles .  Notable  differences 
are:  (1)  Indeterminism  in  classical  statistical  mechanics  merely  reflects 
our  ignorance  of  initial  conditions.  In  the  quantum  world  indeterminism 
seems  to  be  unavoidable  even  in  principle.  (2)  Unlike  classical 
mechanics,  quantum  mechanics  makes  explicit  conflict  between  the  dynamical 
description  and  the  process  of  measurement.  Dynamical  equations  for 
microscopic  systems  are  based  on  objective  reality  and  are  reversible  in 
time,  as  are  the  equations  of  classical  dynamics.  However,  the  observed 
reality  --  the  irreversibility  of  the  microworld  --  is  an  outcome  of 
measurement  and  thus  directly  conflicts  with  the  objective  reality.  (3) 
Bose-EInstein  (BE)  and  Fermi-Dirac  (FD)  statistics  incorporate  quantum 
'tncertainty  through  the  indistinguishability  of  the  particles  with  respect 
to  their  coordinates.  This  gives  rise  to  exchange  terms  in  quantum 
statistical  mechanics.  There  is  no  analogue  of  exchange  terms  in 
classical  statistical  mechanics.  (4)  Many  particle  wavefunctions  for  the 
microworld  obey  certain  symmetry  rules.  This  quantum  mechanical 
limitation  has  no  immediate  counterpart  in  classical  statistics. 

One  possible  way  to  resolve  these  conflicts  would  be  to  abandon  the 
objective  reality,  as  advocated  by  the  believers  of  the  Copenhagen 
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interpretation  [1],  The  quantum  Liouville  equation  for  mixtures,  which 
describes  the  dynamics  of  observed  reality  for  the  microworld,  would  then 
be  an  appropriate  dynamical  equation.  The  central  task  is  thus  the 
evaluation  of  the  quantum  distribution  function  (QDF) ,  because  the  QDFs 
provide  a  means  of  expressing  quantum  mechanical  averages  in  a  form  which 
is  very  similar  to  that  for  classical  averages  [2].  This  is  the  case  with 
the  Ehrenfest  representation  (ER) ,  where  a  correspondence  between  a 
classical  trajectory  and  the  expectation  values  of  the  corresponding 
quantum  operators  is  established  [3],  However,  in  most  cases,  including 
ER,  the  equations  of  motion  for  the  average  values  cf  position,  momentum, 
etc.  are  not  closed  and  therefore  cannot  be  solved  without  further 
assumptions.  Moreover,  attempts  to  write  a  proper  QDF  that  is  positive 
and  gives  the  correct  individual  quantum  distributions  of  position  and 
momentum  have  been  in  vain  [2,4].  It  is  suspected  that  the  quantum 
mechanical  uncertainty  and  the  noncommutativity  of  position  and  momentum 
prohibits  a  phase  space  construction  of  the  QDFs.  For  these  reasons,  it 
has  not  been  possible  to  develop  theories  for  the  microworld  based  on 
observed  reality.  However,  in  recent  years,  it  has  been  found  that  the 
Wigner  transforms  of  all  the  Gaussian  QDFs  are  non-negative,  corresponding 
to  pure  states  and  mixtures  [5].  Existence  of  such  phase  space  QDFs  gives 
us  an  opportunity  to  study  the  dynamics  of  the  microscopic  world  based  on 
observed  reality.  This  is  what  constitutes  the  main  content  of  the 
present  paper. 

Our  system  of  investigation  is  an  N-particle  statistical  system 
whose  time  evolution  is  described  by  the  quantum  Liouville  equation  for 
mixtures  [6].  We  focus  on  studying  the  dynamics  of  each  single  particle 
in  the  N-particle  ensemble.  We  do  not  assume  &  priori  that  the 
indistinguishability  of  the  particles  is  necessary  when  their  de  Broglie 
wavelengths  (A)  are  greater  than  or  equal  to  their  mean  distances  (R) . 
Instead,  we  measure  the  quantum  uncertainties  associated  with  each 
particle  in  order  to  understand  whether  indistinguishability  is  inherent 
in  microscopic  many-particle  systems.  The  development  is  based  on 
constructing  a  Gaussian  density  matrix  in  three-dimensional  phase  space. 
The  density  matrix  is  defined  by  a  set  of  dynamical  variables  whose 
expectation  values  are  considered  to  be  relevant  for  the  dynamics. 
Construction  of  the  density  function  is  based  on  a  .  maximum  entropy 
formulation  [7],  and  our  choice  of  the  dynamical  variables  is  reflected  in 
such  construction.  The  self-consistent  equations  of  motion  are  then 
derived  for  these  expectations  from  the  quantum  Liouville  equation  using  a 
projection  scheme  [8,9].  The  solution  of  these  self-consistent  field 
(SCF)  equations  provides  the  time  evolution  of  the  density  matrix. 

In  the  next  section  we  derive  the  equations  of  motion  for  the 
expectations  of  these  dynamical  quantities  and  construct  the  corresponding 
density  function  for  mixed  states.  We  also  show  that  these  equations  of 
motion  can  be  derived  solving  the  classical  Liouville  equation  for 
mixtures.  In  Sec.  Ill  we  show  that  for  systems  in  thermal  equilibrium 
these  equations  satisfy  a  new  principle  of  least  action.  Again,  a  simple 
combination  of  classical  and  quantum  laws  can  be  used  to  generate  these 
equations.  Finally,  in  Sec.  IV  we  discuss  the  relative  merits  of  the 
theories  based  on  observed  reality  with  those  of  objective  reality. 


DERIVATION  OF  THE  BASIC  EQUATIONS  SOLVING  QUANTUM  LIOUVILLE  EQUATION  AND 
CONSTRUCTION  OF  THE  DENSITY  FUNCTION 


We  characterize  our  N-particle  system  by  the  Hamiltonian 
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1.1  4. 


and  a  density  matrix  p(R,R';t)  which  satisfies  the  quantum  Liouville 
equation  (QLE) 

&  -  •  £[",/>]  S  -iU  ,  (2) 

where  V  is  the  interaction  potential  and  m,  is  the  mass  of„  the  k-th 
particle,  r^  is  the  coordinate  of  the  k-th  particle,  and  p,  -  ~  ^  is  its 
conjugate  momentum.  We  are  interested  in  an  approximate  solution  of  the 
QLE  where  the  time  evolution  of  the  exact  N-particle  system  is  described 
by  a  reduced  density  matrix  /?re(R,R' ;  t) ,  which  is  a  product  of  N  single 
particle  density  functions , 

N 

P  (R,R';t)-n  *.(r.,r;;t)  ,  (3) 

re  i_l  J  J  J 

where  R  is  a  vector  with  N  coordinate  components  r....r  .  This  choice 
[Eq.  (3)]  of  the  density  function  excludes  the  possibility  of 
incorporating  quantum  uncertainty  into  our  system  through  the 
indistinguishability  of  the  particles  with  respect  to  their  coordinates. 
We  introduce  the  effect  of  quantum  uncertainty  by  characterizing  each 
particle  with  respect  to  its  position,  momentum  and  their  fluctuations. 
For  this  reason,  we  define  each  particle  in  three  dimensions  using  a  set 

of  operators,  0  ,  a  -  1 .  M,  which  are  x,,  p,  ,  x,x,,  p.p.  and  x,p.,  for 

k,i  -  1,  2,  3.  In  general,  to  describe  tne  dynamics  of  each  particle  in 

thrtje  dimensions  (3D)  incorporating  quantum  uncertainties,  we  need  M  - 
2.3  +  3.3  -  27  independent  operators  consisting  of  a  complete  set  of 

linear  and  bilinear  products  of  x^  and  p^.  For  notational  convenience,  we 
shall  represent  these  operators  in  the  following  matrix  form: 


-  1  ,  Ax  -  r  ,  A2  -  p 

-  (r  -  ^‘(r  -  ) T  ,  A^ 

-  (i  -  0'1) •  (p  -  a2)T  -  \  MI 


-  (p  -  0*(P  -  Or))' 


Here, 
A.  an 


Aq  is  the  unit  operator;  A.  and  A2  are  three  component  vectors;  A,, 
Ac  are  3x3  matrices  and  I  Is  the  3x3  unit  matrix;  a,  and  are  tne 


A^  and  A^  are  3x3  matrices  and  I  Is  the  3x3  unit  matrix;  a ^  and  are  the 
expectation  values  of  the  dynamical  quantities  A^  and  A2;  and  the 
superscript  T  stands  for  transpose.  The  time  evolution  of  the 
expectations  of  these  dynamical  quantities,  <A.>’s,  will  provide  us  with 
the  average  values  of  position  (<A->),  momentum  (<A~>) ,  and  their 
fluctuations  (diagonal  terms  of  <A~>  and  <A,>)  and  correlations  (off- 
diagonal  terms  of  <A->  and  <A^>  and  <A^>) .  <A^>  nas  the  unit  of  action, 
and  the  role  of  its  diagonal  elements  In  the  present  dynamical  description 
will  be  discussed  later.  Here,  we  should  note  that  <A^>  and  <A^>  are 
symmetric  matrices. 


To  derive  explicit  expressions  for  the  time  evolution  of  the 
expectations,  <0,  (t)>,  for  each  particle  j,  we  use  a  projection  operator 
technique  in  Liouv?lle  space  [8,9].  In  this  space,  ft  and  p(t)  can  be  Codes 
written  as  |H»  and  |p(t)».  Since  the  exact  density  matrix  p(t )  contains  vor 
information  more  than  required  by  the  present  method,  we  evaluate  ^ 
expectation  values  of  the  dynamical  quantities  0ja  using 
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a.  (t)  S  «0.  |p(t)»  -  «0 .  Ip  (t)» 
ja  ja'  ja‘  re 


a  -  0 ,  1 ,  M  .  (5) 


For  each  degree  of  freedom  j ,  we  now  define  an  (M+l)x(M+l)  overlap  matrix 
with  elements 

^  V> 

-  Tr  [0ja  pre(t)  0^]  a,  f)  —  0,  1 .  M  (6) 

and  the  Liouville  space  projection  operators 
N  M 

P(t>  -}  }  «°j*l  •  '7> 

j-1  <.,0-0 

The  properties  of  the  above  projection  operator  have  been  discussed  in 
Refs.  8  and  9.  Let  us  now  assume  that  at  some  time  t  -  t' 

P(t')  -  pre(t')  .  (8) 

Introducing  the  complementary  projection  Q(t)  -  1  -  P(t),  along  with  the 
projection  operator  properties  discussed  in  Refs.  8  and  9,  we  can  then 
derive  the  equations  of  motion  for  the  0.  (t)'s  from  the  QLE  in  the  form 
[9]  i* 

aja(t)  -  -i«OjQ|L|pre<t)»  +  ^  M^(t.t')  a^p(t)  ,  (9a) 

P 

where  M^(t,t')  is  expressed  in  terms  of  the  28x28  matrices 


and  is  given  by 


-  -i«0ja|LQ(t)U(t,t')|pre(t')0.^»  , 

(9b) 

(9c) 

M 

Wi7(t,t')[Rj(t,t')]^  . 

7-0 

(9d) 

the  time  evolution  operator 

exp[ -iL(t-t’ ) ] 

(9e) 

Equations  (9)  describe  the  time  evolution  of  the  27N  dynamical  quantities 

(j  -  1,  2 .  N;  a  -  1,  2,...,  27)  and  are  exact.  These  equations  are 

closed  for  27N  unknown  a,  (t)'s.  In  these  equations,  the  time  derivative 

of  a.  (t)  at  time  t  depends  on  all  a.„(t)  at  the  same  time.  Note  that  we 

assume  0.Q  to  be  the  unit  operator  and  the  normalization  condition 

requires-^  its  expectation  values  to  be  independent  of  time,  a...  -  1.  Now 
if  we  assume  that  condition  (8)  holds  for  all  times,  then  Q(t)p(t)  -  0  and 
we  are  left  with  the  first  term  on  the  right-hand  side  of  Eq.  (9a),  which 
represents  a  mean-field  time  evolution  of  the  N-particle  system,  where  the 
time  evolution  of  the  expectations  of  the  dynamical  quantities,  0.  ,  are 
given  by  J 

4j(j(t)  -  -  J  Tr{OjQ[H,pre(t)]J  E  -  i«Oja|L|pre(t)»  .  (10) 
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Equation  (10)  describes  the  time  evolution  of  the  dynamical  quantities  x,  , 
J>k>  p^  and  xfcpi  for  k,i  -  1,  2,  3.  K 

To  determine  directly  the  quantum  uncertainties  associated  with  each 
particle,  we  evaluate  explicitly  the  time  evolution  of  the  dynamical 
quantities 


a i  -  <r>  ,  a ^  -  <p> 

o3  -  <(r  -  o1)»(r  -  £71)T> 

ff4  "  <(0  *  *  a2)T> 

a ,  -  <(r  -  o,)-(p  -  o„)T>  -  kl 


given  by 


o2  -  -  V' 

*  a  a  T 

a ^  A  +  A 

a.  -  -  V”a,  -  crjv” 
4  5  5 

-  -  c^V"  + 


(11a, b) 
(11c) 
(lid) 
(He) 

(12a) 

(12b) 

(12c) 

(12d) 

(12e) 


where 


dr-dr 


A  -  m  ^1 


Equations  (12)  describe  time  evolution  of  27  matrix  elements  which  are  the 
expectations  of  the  physical  operators  described  by  Eq.  (4).  a  for  a  - 
1 ,  2 .....  M  are 


fflk(t)  -  <V  ■  a2k(t)  “  <Pk> 

(14a, b) 

[cr3(t)  Jki 

-  <[\  •  ‘,ik<t)llx»  ■ 

*„«:>]> 

(14c) 

[^4(t)]ki 

'  <[pk  ‘  <’2k(t)i!l’i  • 

(14d) 

U1 
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s/ 

Ka 

"  2  <i*k  ■  ,’lk<t)IIpl 

+  [p«  '  ’2i(t)llSk  ■  1 

’lk(t)1> 

(14e) 

k,i  -  1,  2,  3 

These  expectation  values  in  quantum  mechanics  are  evaluated  either  in  the 
coordinate  (r,r')  or  momentum  (p,p')  representation. 

From  Eq.  (12),  we  find  that  for  a  successful  application  of  the 
present  method,  the  choice  of  4>  is  crucial.  We  construct  each  d,  from 
the  physical  consideration  which  ii  the  maximum  entropy  principle  [ 7-j  .  We 
associate  entropy  with  each  particle  using  the  relation 
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We  then  maximize  entropy  (5S  -  0)  subject  to  the  constraint  condition 


Saa(t)  -  Tr[0a  6^]  -  0 


a  -  0,  1 .  M  . 


(15b) 


This  yields  the  maximum  entropy  distribution 


*J<t>  "  6XP  [  I  '  Aja(t)0jJ  • 


where  the  X.  (t)'s  are  Lagrange  multipliers  and  describe  the  time 
evolution  of  the  density  function  <j>  .  This  distribution  function  is 
Gaussian  in  30-phase  space  (q,p)  ^and  in  the  corresponding  coordinate 
representation  (r,r'),  whose  explicit  forms  are  given  by 


<£(r,r';t)  -  C  exp{-  ^(^(r+r')T  o^Cr+r')  -  rn  (r-r '  )T(<7^a^1a1-a,) 


+  ^2  (r-r' )T(o4-ff^or‘1o5)  (r-r' ) 
Y- 


+  i^  ^r+r')T  a31<75(r*r')  -  (r+r' )Ta31o1] ) 


(17a) 


(2*f) "3//2  ||o3|C1/2  exp( -  ^■1a1I  , 


(17b) 


tf(q,p;t)  -  ( 2n )  3  ||B||*I/2  exp{-  (qT-a*)  ,  (p1-^)  ]  B-1^  .  )  . 

(?8a) 


where  B  is  a  matrix  given  as 


■-r?  t  • 

/T  -* 


(18b) 


^(q,p;t)  is  obtained  form  ^(r,r';t)  using  the  Wigner  transform  [2] 


*(q.p;t)  -  [[[  ds  <q-s|$(t) |q+s>e2ip’S^ 


(19a) 


and  satisfies 


JJ<*’ 


dp  ^(q,p;t)  -  1 


(19b) 


Therefore,  the  expectation  values  of  the  dynamical  quantities  [Eq.  (11)] 
in  the  present  method  can  be  evaluated  using  both  the  trace  operation 
(r,r')  and  the  phase  space  (q,p)  integration  methods.  The  time  evolution 
of  these  density  functions  are  obtained  by  solving  the  coupled  set  of 
equations  (12).  A  close  look  at  these  equations  shows  that  they  do  not 
contain  Y-  This  suggests  that  they  are  completely  classical.  In  fact,  we 
obtained  these  TDSCF  set  of  equations  (12)  also  by  solving  the  classical 
Liouville  equation  using  the  phase  space  distribution  function  (18),  and 
following  the  same  procedure  discussed  above.  The  TDSCF  set  of  equations 
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The  present  classical  set  of  equations  (12)  differ  from  Ehrenfest's 
representation  due  to  the  fact  that  they  contain  additional  dynamical 
quantities  (o^,  a ^  and  and  are  closed.  To  understand  the  potentiality 
of  the  present  method,  we  need  extensive  application.  So  far,  we  have 
studied  only  the  equilibrium  properties  of  liquid  and  solid  argon, 
assuming  the  system  to  be  isotropic.  The  results  of  this  simulation  have 
been  reported  in  a  separate  communication  [10].  Here  we  elucidate  only 
the  underlying  physics.  Of  the  various  dynamical  quantities,  only  the 
time  evolution  of  orj(t),  which  has  the  unit  of  action,  shows  rapid 
oscillations  between  positive 


and  negative  values.  This 
time  dependence  of  a^( t) 
implies  that  at  different 
times  of  the  measurement  each 
particle  has  a  different 
trajectory,  which  in  turn 
suggests  that  the  classical 
trajectory  obtained  solving 
(12)  does  not  satisfy  the 
principle  of  least  action 
[11].  The  time  evolution  of 
this  (7e ( t)  is  shown  in  Fig.  1. 
From  tnis  figure,  we  find  that 
the  time -averaged  trajectory 
for  systems  in  thermal 
equilibrium  has  -  0.  Since 
our  system  of  investigation  is 

i  *  j 

times . 
in 


isotropic,  [ <7^5 j  ^ . 
remains  zero  ovei~ 
Now  substituting 
Eq.  (12), 


for 
all 
<A5>  -  0 
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over  20  p s  in  a 


of 
for 
typical 


1  -  x. 


action, 

v  and  z 

run  for 
50.275 


solid  argon  (V  -  28.096,  .  -  - - 

K) .  ffi( t)  is  the  action  along  the 
i-th  direction  of  phase  space  and 
averaged  over  the  entire  bath.  The 
value,  of  o-(t)  originally^, in  kg 
m"  s  has  been  magnified  10  times 
in  this  figure. 


That  is,  for  systems  in  thermal  equilibrium,  the  fluctuations  and 
are  the  constants  of  motion.  In  the  following  section,  we  show  tnat  Eq. 
(20)  can  be  obtained  using  a  modified  representation  of  the  classical 
principle  of  least  action  and  also  from  a  simple  combination  of  the 
Ehrenfest  representation  with  the  equations  of  motion  derived  from  the 
quantum  uncertainty  principle. 


ALTERNATIVE  APPROACHES  TO  THE  DERIVATION  OF  THE  BASIC  EQUATIONS  FOR 
SYSTEMS  IN  THERMAL  EQUILIBRIUM 

Principle  of  Least  Action 

We  introduce  the  following  principle  of  least  action  for  each 
particle  in  the  N-particle  ensemble: 


8 


-  Jt 


% 


dt  <Lj> 


j  -  1,  2 .  N 


(21) 


Here  is  the  Lagrange  operator  for  the  j-th  particle  and 


3  a2 


j  Z  2m  +  Vj 

k-1  -> 


(R) 


(22) 


where  its  quantum  average  is 

3 

Lj  -  <V*  2±  ]  l(°2k>2  *  KWj  -  f  <avj(R)  Pre(R,R;t)  . 

J  k-1  *® 

(23) 

and  the  o^'s  are  defined  in  Eq.  (14).  Let  us  ipw  assume  that  the 
extremum  path  f^r  the  i-th  particle  is  defined  by  cr*(t),  whose  end  points 
are  °r^(ta)  and  ^(t^)-  In  3D  phase  space 


9  2  l  i  i  i  1 

-  i<x  >  +  j<y  >  +  k<z  >  -  io^  +  jo^  +  ^ 


13 


calculus  of  variations . 


-2 

£  Let  us  assume  that  the  path  is  varied  away  from  cr.(t)  by  an  amount 
So^t),  with  the  condition  that  the  end  points  are  fixed,  that  is, 


«ff{(ta)  -  ^(V  -  0  . 

-2 

Now  the  condition  that  o^(t)  be  an  extremum  means 

,  l  l.-i  ^  c  i,  i .-i.  . 

oo^  -  +  o*7^]  "  "  0 


(25) 


(26) 


to  first  order  in  So Using  the  definitions  (21)  and  (23),  we  may  write 

t, 

ir  i  ^  ,  i,  f  r  /  *  j.  r  *  i  f  i  i  ^  -  i  i  ^  .  2. 

ff5^al  +  SaV  ~  J  dt  Li^l  +  5al-  a2  +  S<J2’  a3  +  &<3Z'  ff4  +  5cV 


i,-. 

“  o,.[o 


2  2  i 

{]  +  J  dt  {So;  f  +  So1  ~ 


1  .  i 
do^ 


2  .  i 
do  2 


2  ^L2  l  3Lf 

+  s°l  Tm  +  5\  TV 

do~  do, 

3  4 


(27) 


where ^  for  example,  dL./dof  implies  variation  of  L£  with  respect  to  all 
the  elements  (k-1,  2,  3;.  Upon  integration  by  parts,  the  variation 

r  D« 


in  becomes 


i  i  3Li 

5o5  -  m£  S ox  i 

do  2 


^  A  .  j  3L. 

-  j  dt  {5o1[mi  ^(-f) 

C  * C  60 ^ 

a  a  2 
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| 
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dL.  .  3L.  .  3L» 

-  — j]  +  So*  — j  +  So*  — J)  .  (28) 

do ^  3£7j  3ff^ 

l 

The  first  term  on  the  right-hand  side  of  this  equation  is  zero  since  So ^ 
is  zero  at  the  end  points.  The  second  term  is  zero  if  the  following 
condition  is  satisfied: 

i  d  dLg  dLg  2  2  ^2 

‘  2]  +  5*3  i  +  6al  2  ~  0  •  (29) 

1  az  do*  a o*  So*  4  3o£ 

i  i  k  k 

Since  the  variation  of  with  respect  to  o ^  and  contains  and  <73  for 
k  *  2,  to  obtain  the  extremum  path  for  each  particle  we  need  to  solve  3N 

coupled  equations  obtained  from  the  variation  of  ff,  (2  -  1,  2,...,  N) .  In 
the  following,  we  shall  show  that  Eq.  (20)  can  result  from  a  trivial 
solution  of  these  3N  coupled  equations.  Since  between  the  end  points  5o^ 
can  take  any  arbitrary  value,  the  first  term  in  Eq.  (29)  is  zero  provided 

d  3Li  dh2 

■  71  “ 0  ■  (30) 

do ^  do^ 

These  are  the  Langr^nge  equations  of  motion  for  the  quantum- averaged 
dynamical  quantities  o.(t)  and  <r_(t).  The  second  and  third  terms  in  Eq. 
(29)  can  only  be  zero  If  * 

$03  -  SoZu  -  0  .  (31) 

2 

T^is  is  because,  is  we  differentiate  of  Eq.  (23)  with  respect  to  <7  3  and 

0 .  ,  then 
4 


—7  -  -  ,  <V„(R)>  *  0  ,  ~  -  .  (32) 

dot  dot  1  do*  m2 

3  3  4 

Solving  the  Lagrange  equations  (30),  we  obtain  Eqs .  (20a, b) ,  and  from  the 
condition  (31)  we  obtain  Eqs.  (20c,d).  Thus,  a  trivial  solution  of  the  3N 
coupled  equations  shows  that  there  is  a  principle  of  least  action  (21) , 
from  which  Eq.  (20)  can  be  obtained.  In  the  following  subsection,  we  show 
that  Eq.  (20)  can  be  derived  by  combining  the  Ehrenfest  representation 
with  the  equations  of  motion  for  the  fluctuations  obtained  from  the 
quantum  uncertainty  principle. 


Ehrenfest' s  equations  describe  the  time  evolution  of  the  mean  values 
of  the  coordinates  <x>  and  momenta  <p>.  With  the  Hamiltonian  given  by  Eq. 
(1),  the  time  evolution  of  these  expectations  may  be  written  as 


^  -  <pj/m>  j  - 1 .  2 , . . . 

fc"j>  -  ■  ^  ■  <vj(R)>  J-1’ 


,2 .  3N  . 


(33a) 


(33b) 


These  equations  do  not  follow  the  classical  laws  of  motion  due  to  the  fact 
that  the  time  derivatives  of  <x,>  and  <p,>  are  equal  to  certain  average 
values  whose  calculation  generally  necessitates  the  knowledge  of  the 
wavefunction  ♦(t).  Our  TDSCF  equations,  (20a, b),  are  identical  to  these 
equations,  except  that  the  average  values  are  evaluated  using  a  density 
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function  p  (t).  Ehrenfest's  equations  may  thus  be  viewed  as  a  special 
case  of  our  TDSCF  set  of  equations  (20) ,  where  the  contributions  of  the 
fluctuations  are  not  taken  into  account  and  the  density  function 
represents  a  pure  state  [3,6], 

Although  the  TDSCF  set  of  equations^ (20)  arj  classical,  2che  timj 
evolution  of  the  fluctuations  a ,  -  <x^>  -  <x.>  and  -  <p.>  -  <p,> 

shows  that  they  always  satisfy  the  ^quantum  uncertainty  condition.  The 
time  evolution  of  a ,  ,  and  a,  .  for  liquid  (V  -  29.283,  T  —  89.834  K)  and 
solid  (V  -  28.096,  x  ^  50.275^K)  argon  usingmthe  Leonard- Jones  potential 
are  shown  in  Figs.  2  and  3,  respectively. 
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Fig.  2. 


Time  evolution  of  the  dynamical  variable  t)  - 


3 

2  c7 «  ( t )  in  the 

-22  2  **  i*»l  ^ 

units  10  m  obtained  by  solving  Eq.  (20)  using  the  predictor 

algorithm  for  liquid  (---)  and  solid  (-)  argon  at  the 

temperatures  89. 834. K  and  50.275  K,  respectively,  and  over  a  time 

period  of  40  ps.  Oj(t)  represents  the  position  fluctuation  along 

the  i-th  direction  of  phase  space  and  averaged  over  the  entire 

bath. 


Fig.  3.  The  same 


caption  as  for  Fig.  2,  but  replacing 

46  2  2  2 

vertical  axis  with  the  units  10  kg  m  s 


<7^  with  a ^ 


on  the 


In  general, 
satisfy 


a3Aj 


we  observe  that  the  magnitudes  of  the  fluctuations  always 


w2 

-Q>-\ 


j  -  1,  2, 


3N 


(34) 


where  the  magnitude  of  G  is  directly  proportional  to  the  energy  (<E>)  of 
the  system.  This  in  turn  suggests  that  if  we  have  a  system  such  that  the 
fluctuations  satisfy  condition  (34),  then  its  time  differentiation  yields 
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3N 


(35) 
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Since  and  a are  not  zero  from  the  uncertainty  condition,  one  of 
possible-'ways  to-1  satisfy  (35)  would  be 


the 


3j 


j  -  1,  2... 


3N 


-  0 


j  -  1,  2,...,  3N 


(36a) 

(36b) 


Combining  Eqs.  (33)  and  (36),  we  obtain  a  set  of  equations  analogous  to 
(20) ,  except  that  the  choice  of  the  density  function  is  still  open  in  the 
present  case.  The  density  function  may  be  for  a  single  particle  or  for  an 
N-particle  system.  Again,  the  density  function  may  represent  a  pure  state 
or  a  mixture.  For  a  Gaussian  choice  of  the  density  function,  our  TDSCF 
set  of  equations  is  thus  capable  of  describing  the  dynamics  of  all 
macroscopic  and  microscopic  systems. 


DISCUSSION 


of  observed 
statistical  , 
in  phase 


is 


We  have  described  a  method  for  studying  the  dynamics 
reality.  The  observed  system  may  be  classical,  classical 
quantum,  or  quantum  statistical.  The  present  development 
space  and  thus  recovers  the  trajectory  concept  even  for  the  microscopic 
systems.  In  the  present  method,  the  trajectory  of  a  particle  is  described 
by  a  coupled  set  of  equations  (12) .  These  equations  represent  the  time 
evolution  of  the  mean  values  of  position  (cr.,)  and  momentum  (o.)  ,  their 
fluctuations  (diagonal  elements  of  a ^  ana  a.) ,  their  correlations  (off- 
diagonal  elements  of  a.,  a ^  and  a^) ,  and  the  action  terms  (diagonal 
elements  of  ct^).  The  present  study  has  been  restricted  to  isotropic 

systems,  where  off-diagonal  elements  vanish. 

The  diagonal  elements  of  <7,.(t)  are  identified  with  action  terms 
since  they  have  units  of  action  and  their  time  evolution  is  similar  to 
that  of  the  phase  terms  (o^/)tf)  in  the  path  integral  approach.  Thus,  our 
present  method  is  reminiscent  of  Feynman's  path  integral  approach  [11], 
where  different  values  of  action  correspond  to  different  trajectories.  In 
the  path  integral  method,  optimal  trajectories  are  obtained  by  invoking 
the  classical  principle  of  least  action.  Inspired  by  this  approach,  we 
have  introduced  a  modified  version  of  the  classical  principle  of  least 
action  [Eq.  (21)]  in  the  analysis  of  our  results.  It  is  worth  noting 
that,  in  the  present  method,  the  density  function  contains  phase 
information.  In  contrast  to  the  path  integral  approach,  the  phase 
vanishes  over  the  minimum  action  path. 


For  a  Gaussian  choice  of  the  density  function  [Eqs.  (17)  and  (18)], 
the  equations  of  motion  (12)  are  closed  and  are  completely  classical.  For 
an  N-particle  statistical  system,  we  derived  these  equations  from  the 
quantum  Liouville  equation  using  a  projection  scheme.  It  was  also  assumed 
that  the  SCF  density  function  approximates  the  exact  density  function  over 
all  time.  This  choice  of  the  density  function  as  given  in  Eq.  (3)  is 
similar  to  the  Hartree  approximation,  where  indistinguisability  of  the 
particles  with  respect  to  their  coordinates  is  ignored. 

The  concept  of  indistinguishability  in  quantum  statistics  did  not 
originiate  from  first  principles,  but  rather  from  thj  need  to  explain 
quantum  uncertainty.  Indeed,  that  indistinguishabilty  is  not  inherent  in 
quantum  statistics  is  evident  in  our  model.  We  characterized  each 
particle  with  respect  to  its  position  and  momentum,  and  their  respective 
fluctuations.  The  time  evolution  of  the  fluctuations  allows  direct 


determination  of  the  quantum  uncertainties  associated  with  each  particle 
(Figs.  2  and  3).  These  fluctuations  are  constants  of  motion  over  the 
minimum  action  path  ( So ^  -  0,  o>  -  0)  and  satisfy  the  quantum  uncertainty 
condition  (34).  This  suggests  that  indistinguishability  is  not  a  must  in 
quantum  statistics.  Thus,  all  that  is  needed  to  describe  an  N-particle 
system  is  N  trajectories  with  the  constraint  of  minimum  action. 

It  may  seem  that  the  validity  of  this  classical  description  (20) 
results  from  the  choice  of  a  normalized  Gaussian  density  function  [Eqs. 
(17)  and  (18)],  and  therefore  the  inclusion  of  higher-order  moments  would 
cause  a  departure  from  the  classical  picture.  By  taking  higher -o^der 
moments,  we  introduce  fluctuation  terms  such  as  a,  -  <x,>  -  <x.Xx;>  - 
<x,>  ,  which  are  of  third  order  or  greater.  Of  course,  these  higher-order 
terms  represent  a  departure  from  the  Gaussian  picture,  but,  at  the  same 
time,  go  beyond  the  variances  (a^ , a^)  defining  uncertainty  in  quantum 
measurement.  Insofar  as  quantum  measurement  is  random,  the  choice  of  a 
Gaussian  density  function  is  appropriate. 

Computer  simulation  shows  that  the  fluctuations  in  position  and 
momentum  of  each  particle  are  constants  of  motion  over  the  minimum  action 
trajectory.  Thus,  it  appears  that  each  particle  maintains  its  distinct 
trajectory  over  all  time  even  though  the  de  Broglie  wavelengths  of 
individual  particles  overlap.  Insofar  as  computer  simulation  suggests 
that  there  are  N  distinct  trajectories  for  an  N-particle  system,  the  SCF 
approximation  (3)  is  vindicated. 

The  validity  of  the  SCF  approximation  in  turn  suggests  that  ai.  N- 
particle  system  is  analogous  to  N  interacting  wave  vectors  in  a  Hilbert 
space.  Measurement  on  each  wavevector  then  yields  the  density  function  <j>. 


<f>  -  )  P  |u  >  <u  I 
/  n  1  n  n1 


From  the  viewpoint  of  quantum  mechanics,  the  measurement  of  a  dynamical 
quantity  gives  various  eigenvalues  with  corresponding  eigenvectors  |u  > 
and  probabilities  P  .  Measurement,  therefore,  introduces  irreversibility 
in  our  microscopic  description.  Indeed,  trajectories  obtained  from  the 
Schrddinger-Heisenberg  picture  are  indeterministic  and  consequently  time- 
irreversible.  Yet,  a  phase  space  (q,p)  representation  of  the  density 
function  in  Gaussian  form  (18)  recovers  the  element  of  time  reversibility 
of  measurement  (20).  Notably,  computation  of  trajectories  in  our  method 
requires  knowledge  not  only  of  the  position  and  momentum  but  also  of  their 
fluctuations . 

Computed  trajectories  in  the  present  method  are  unique  since  the 
density  function  chosen  is  the  one  satisfying  the  maximum  entropy 
principle.  This  distribution  has  an  important  advantage  over  the 
canonical  distribution  (e '^H)  in  that  the  energy  constraint  condition 
(i<H>  -  0)  is  eliminated.  Our  choice  of  dynamical  constraints  (15b)  in 
the  construction  of  the  density  function  (15a)  makes  our  theory  suitable 
for  nonequilbrium  statistical  processes.  Thus,  the  present  method  is  far 
more  general  than  the  equilibrium  formulations. 

In  the  classical  ensemble  picture,  individual  configurations  have  no 
significance,  whereas  in  the  quantum  case  it  is  possible  to  obtain  energy 
eigenvalues  (instead  of  average  energies)  for  bound  state  systems,  even 
though  the  time-energy  uncertainty  holds.  This  is  because  the  uncertainty 
principle  imposes  restriction  only  on  the  product  of  the  uncertainties  in 
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position  and  momentum  (or  equivalently  energy  and  time)  measurement  and 
not  on  the  individual  uncertainties .  Relying  on  this  premise  and  noting 
that  bound  states,  because  of  their  standing  wave  nature,  have  no  explicit 
time  dependence,  the  existence  of  individual  bound  energy  eigenstates  is 
assumed  in  the  Schrddinger  picture.  Solution  of  the  corresponding  wave 
equation  gives  the  energy  eigenvalues.  A  striking  departure  from  the 
classical  ensemble  picture  is  thus  evident  in  quantum  physics.  In  the 
present  development  the  concept  of  the  ensemble  picture  is  retained  even 
in  the  construction  of  the  density  function.  Therefore,  the  concept  of 
eigenstates,  which  is  central  to  BE  and  FD  statistics,  is  superfluous  in 
the  present  method.  Further,  the  symmetry  rules  of  Pauli  are  implicit  in 
our  approach  since  we  do  not  need  explicit  knowledge  of  the  probabilities 
P  [Eq.  (37)].  In  order  to  adapt  the  present  method  to  Fermions, 
additional  variables  are  necessary  to  describe  the  spin  orientation 
associated  with  each  particle.  This  extension  will  be  reported  in  the 
future . 

So  far,  we  have  confined  applications  of  the  present  method  to 
statistical  systems  only.  Nevertheless,  it  is  also  suitable  for  studying 
atomic  and  molecular  systems.  S^nce  for  an  N-particle  system  we  need  a 
priori  knowledge  of  N  rather  than  N  density  functions,  our  approach  is 
particularly  attractive  for  large  molecular  systems.  These  applications 
await  future  studies.  In  the  present  method,  the  kinetic  energy  (KE)  is 
evaluated  from 
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The  temperature  for  a  classical  system  is  related  to  the  kinetic  energy  by 
the  equipartition  theorem, 
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Because  of  the  appearance  of  the  momentum  fluctuation  terms  a the 
appropriate  definition  of  temperature  in  our  approach  is  yet  to  be 
resolved.  Here  we  should  note  that  temperature  in  Figs.  2  and  3 
incorporates  the  effect  of  fluctuation. 

In  summary,  the  present  dynamical  description  requires  knowledge  of 
position,  momentum  and  their  respective  fluctuations.  This  does  not  allow 
precise  measurement  of  either  position  or  momentum.  As  a  consequence,  the 
Einstein,  Podolsky  and  Rosen  paradox  [1,4]  does  not  arise. 
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